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Abstract 

We investigate tunneling properties of Bogoliubov mode in a Bose-Einstein condensate. Us- 
ing an exactly solvable model with a ^-functional barrier, we show that each component in the 
two-component wavefunction (n, v) of low-energy Bogoliubov phonon has the same form as the 
condensate wavefunction in the supercurrent state. As a result, the currents Ju and associated 
with u and v, respectively, have the same tunneling properties as those of supercurrent carried by 
condensate. Thus, the tunneling of low-energy Bogoliubov phonon described by the tunneling of 
these two currents shows perfect transmission. We also show that the supercurrent behaviors of 
Bogoliubov phonon still exist in the presence of supercurrent carried by condensate, except in the 
critical supercurrent state. In the critical current state, the perfect transmission is absent, because 
Ju or Jy exceeds their upper limit given by the critical value of the supercurrent associated with 
the condensate. Our results consistently explain the recently proposed two tunneling phenomena 
associated with Bogoliubov phonon, namely, the anomalous tunneling effect (perfect transmission 
in the low-energy limit) and the breakdown of the perfect transmission in the critical supercurrent 
state. 
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I. INTRODUCTION 



Recently, novel tunneling phenomena of Bogoliubov phonon have been theoretically pre- 
dicted in superfluid Bose gases. Kovrizhin and co-workers [l, 2, sl clarified the perfect trans- 
mission of low-energy Bogoliubov phonon across a potential barrier, which is referred to 
as the anomalous tunneling effect. Danshita and co-workers^] showed that the anomalous 
tunneling effect also occurs in the supercurrent state, as far as the magnitude of the su- 
percurrent is less than the critical current |^]. In the critical supercurrent state, the perfect 
transmission is not obtainedj4|, irrespective of the relative direction between the momentum 
of Bogoliubov phonon and superflow. They also extended their work to a Bose condensate 
in an optical lattice [s], Q]. Although the anomalous tunneling effect has not been observed 
yet, a cold atom gas may be useful to examine this interesting tunneling phenomenon. With 
this regard, we briefly note that a double- well trap has been realized in cold atom physics [3], 
and a kind of Josephson effect has been observed js]. 

For the mechanism of the anomalous tunneling effect, various key issues have been dis- 
cussed, such as resonance tunneling[3], localized state near the potential barrier [4 1, similarity 
between the wavefunction of Bogoliubov mode and condensate wavefunction in the low en- 
ergy limits], and coupling of quasiparticle current with supercurrent near the barrier [lo|. 
However, despite these great efforts, no consistent explanation for the anomalous tunneling 
and the breakdown of this effect in the critical supercurrent state has not been given yet. 
Since the appearance of Bogoliubov phonon is one of the most fundamental phenomena 
in the superfluid phase [l^, clarifying physical properties of this collective mode is a very 
important issue in the research of superfluidity. 

In this paper, we investigate tunneling properties of Bogoliubov excitations in a weakly in- 
teracting Bose superfluid at T = 0. We treat the condensate wavefunction and Bogoliubov 
excitations within the Gross Pitaevskii (GP) equation and Bogoliubov equations, respec- 
tively. Using an exactly solvable model, we analytically show that the tunneling mechanism 
of low-energy Bogoliubov phonon is the same as that of ordinary supercurrent associated with 
the condensate. Since the supercurrent is well known to tunnel through a barrier without 
reflection, our result immediately explains the anomalous tunneling effect associated with 
Bogoliubov phonon. However, in contrast to the ordinary supercurrent, the Bogoliubov 
phonon consists of two current components (= Ju and J„), whose directions are opposite to 
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each other. These counterpropagating currents have been recently observed experimentally 



by using Bragg spectroscopy 111]. In this paper, we show that they have the same upper limit 
which equals the upper limit of the ordinary supercurrent (critical supercurrent Jc). In the 
critical supercurrent state, J„ or J^, is shown to always exceed J^? so that the supercurrent 
behavior of Bogoliubov mode (perfect transmission) is destroyed. This result gives simple 
and physical explanation for the breakdown of the anomalous tunneling effect in the critical 
supercurrent state predicted in Ref. 

This paper is organized as follows. In Sec. II, we present our tunneling model with a 
5-functional potential barrier. In this model, exact solutions for the GP and Bogoliubov 
equations have been derived in Refs. Since we use these solutions in later sections, 



we summarize them in this section. For their detailed derivations, we refer to Refs. 3, 4|. 
In Sec. Ill, we consider the case in the absence of the supercurrent. Here, we show that 
the anomalous tunneling effect can be explained as a result of the supercurrent behavior of 
Bogoliubov mode. This result is extended to the supercurrent state in Sec. IV. In Sec. V, 
we discuss how the breakdown of the anomalous tunneling effect in the critical supercurrent 
state can be understood based on our results obtained in Sees. Ill and IV. Throughout this 
paper, we set /?. = 1. 

II. MODEL TUNNELING PROBLEM, CONDENSATE WAVEFUNCTION AND 
BOGOLIUBOV EXCITATIONS 

We consider a superfiuid Bose gas at T = 0, which is separated by the potential barrier 
V{x) = V6{x) (Fig. [1]). Since the system is uniform in the y- and z-direction, our model is 
essentially a one-dimensional system. Thus, we only retain the x-direction, and ignore the 
y- and z-direction. We also ignore effects of a harmonic trap potential, for simplicity. The 
latter simplification is allowed when a box-shaped trap is considered [13|. 

In this model, exact solutions for the GP equation and Bogoliubov equations have been 
derived in Refs. 12, 41. Since we use these solutions in later sections, we summarize them in 
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this section. For detailed derivations of the exact solutions, we refer to Refs. [21, IJ| . 



The Gross-Pitaevskii (GP) equation for the condensate wavefunction \E'(x) is given by IJ] 
(-^^ - /i + V6{x) + gmx)\')nx) = 0. (2.1) 
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FIG. 1: Schematic picture of our model. A Bose superfluid is separated by the barrier potential 
V5{x), and we consider the tunneling of Bogoliubov phonon injected from x = — oo. We examine 
both the cases with and without supercurrent Jg (> 0) carried by condensate. The quasiparticle 
current Jqp carried by Bogoliubov phonon is described by the sum of two currents as Jqp = Ju + Jv, 
while the probability current density J is given by J = — J^, where J„ and Jy are defined by 
Eq. (j2.23p . We note that the flow direction of J„ is opposite to that of Ju- 

Here, m is the mass of a Bose atom, and fi is the chemical potential, g is a. repulsive 
interaction between Bose atoms. Introducing the scaled variables, 

^(x) u - V 



V = 



we can write Eq. (12. ip in the dimensionless form 

1 



X 



(2.2) 



2dx^ 



-Ji + VSix) + |^(x)n^(x) = 



(2.3) 



In Eq. (12. 2p . no = |\l/(x = ±oo)p is the condensate density far away from the barrier, 
and ^ = 1/ y/mgno is the healing length. The scaled chemical potential jl equals unity 
in the absence of supercurrent. In the supercurrent state with momentum q, one finds 
/i = 1 + g^/2[4]. In the following, we simply write fi, V, x) as (\E', /x, V, x). 
The boundary conditions at a; = are given by 



^(+0) = ^(-0), 



(2.4) 



dx J x=+o \dx J x=-o 



2r^(0). 



(2.5) 
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We solve the GP equation (12.31) together with the boundary conditions in Eqs. (12.40 and 
(12. 5p . The condensate wavefunction '^q{x) in the supercurrent state is obtained as|4| 



J[qx-sgn(x) 



'''[7(0;) — igsgn(x)], 



where 



7(x) = a/1 - g2tanh(A/l - q'^[\x\ + Xq]), 
ie„ _ 7(0) - iq 



xo in Eq. (12.70 is determined by 



7(0f + V (7(0)2 + g^) - (1 - g2)7(0) = 0. 



(2.6) 

(2.7) 
(2.8) 

(2.9) 



Equation (12.91) is obtained from the boundary condition in Eq. (12.51) . 

We note that the magnitude of the condensate wavefunction \I/g(x) in Eq. (12.61) is sup- 
pressed near the barrier. However, the supercurrent density Js is independent of x as, 



Jo = Im 



(2.10) 



Namely, the supercurrent is conserved in the whole system. 



For a given condensate wavefunction ^I/g(a;), the two-component wavefunction {u, v) of the 
Bogoliubov mode is obtained from the Bogoliubov equations[l^. Using the scaled variables 
in Eq. (12. 2p . one can write the Bogoliubov equations in the dimensionless forms 



'2d^ - + V^i^) + 2|*g(x)|2jM(x) - ^^gixfvix) = Eu{x), 

1 



(2.11) 



{~2d^ - + V6ix) + 2\^g{x)\^y{x) - mxfu{x) = -Eu{x). (2.12) 

Here, the Bogoliubov wavefunction (u{x),v{x)) and the energy E have been scaled as 
{U{x),v{x)) = {u{x),v{x)) / ^/riQ and E = E/gno, respectively. In the following, we omit 
the bars in the Bogoliubov equations (12.111) and (12.121) . The boundary conditions for the 
Bogoliubov mode {u, v) at x = are given by 

'du{x)\ fdu{x) 



«(+0) =m(-0), 
t;(+0) = t;(-0). 



dx / x=+a 
dv{x) 



dx ) x=+o 



dx / x=-o 
dv{x) 



dx / x=-o 



2Vv{Q). 



(2.13) 
(2.14) 
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To construct {u,v) satisfying the boundary conditions in Eqs. fl2.13p and fl2.14p . we need 
particular solutions of the Bogohubov equations (12. lip and (12.120 for x > and x < 0. 



For a given condensate wavefunction ^^^(a;) and energy E, the coupled equations (12.111) and 



(I2.12P have four particular solutions {un,v. 

^2 



Ur 



^i[{pn+q)x-sgn{x)eq] 



1 + ^ ) 7(a^) - isgn{x) { q + 



2J/ 



{n = 1,2,3,4), given by 

Pn I 



a eq u 

0,a 



2E 



l-||)7(x)+^sgn(x)(g + ||(l-g^ 



(We note that Eq. (I2.15P is not normalized.) The momenta p„ {n 
from the expression for the Bogoliubov excitation spectrum|4] 



^ = Pn?+\/Y(Y + 2). 



(2.15) 



1, 2, 3, 4) are obtained 



(2.16) 



Among the four solutions, two of them (n = 1 and 2) describe propagating wave charac- 
terized by real momenta {pi and p2)- The remaining two solutions (n = 3 and 4) describe 
localized states having complex momenta (ps and pi). While only the propagating solutions 
are necessary in considering a uniform system, one has to also take into account the local- 
ized solutions in the present inhomogeneous system. Indeed, in Sec. Ill, we show that the 
localized states appear near the barrier. 

In the low energy region {E <^ 1), p„ {n = 1,2,3,4) reduce tojl^, within the accuracy 
of 0(E), 

E 



Pi 

P2 
P3 
P4 



1 + 



1-q' 



2t^l-q^ + 



-1i^\-q^^ 



qE 
qE 



(2.17) 



1 

Using these particular solutions, we construct the Bogoliubov wavefunction (u,w). As- 
suming that the incident Bogoliubov phonon comes from x = — oo, we set 





-X + 



M>(x) 

t;>(x) 



^(x). 



(2.18) 
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Here, 9{x) is the step function, and 




"^^^^ \ +A \ 1+5 1 ^^^^^ 



{x < 0), 



(2.19) 





+ D I ""'^"^^ I (x > 0). (2.20) 

V3{x) 



The coefficients (A, B, C, D) are determined so that the boundary conditions in Eqs. (12.13^ 
and f l2.14p can be satisfied. We will give their detailed expressions in later sections. 

Once (y4, B, C, D) are determined, the transmission probability is obtained from the ratio 
of the probability current density J = Ju — Jv for the incident wave to that for the transmitted 
wave. Here, J„ and Jy are given by [if 



J u 

J V 



Im 
— Im 



uixY—uix) 
d 

vix)* —vix) 
ax 



(2.21) 



We note that the quasiparticle current density Jqp carried by Bogoliubov phonon has the 
different form from the probabihty current density as Jqp = Ju + Jv To see the difference 
between J and Jqp, it is useful to consider a uniform system, in which the Bogoliubov 
equations give the plane wave solution [l^ 



Ur. 



1 
2 



E 



+ 1 e 



1 + 1 



(2.22) 



E 



1 e 



Here, E = {p^ / 2)[p^ / 2 + 2) is the Bogoliubov excitation spectrum. Substituting Eq. 
f l2:22|) into Eq. ([Mil), we obtain (p > 0), 

'"^'I'^'^i). (>0), ^ ^ 

' (2.23) 

i)p (<o). 



Ju 



J 7J 



2 V E 
1 fp^/2 + 1 



2 V E 

Equation (12.231) shows that the leading terms of J„ and with respect to p are constant in 
the low energy limit. While they dominantly contribute to the probability current density 
J, they are cancelled out in the quasiparticle current Jqp, because the flow directions of Ju 
and Jy are opposite to each other. The contribution to Jqp comes from higher order terms 



in Ju and in terms of p. In Sec. Ill, we find that the cancellation of the leading terms of 
Ju and also occur in the presence of the barrier, as schematically shown in Fig. [H Thus, 
in considering the tunneling of low-energy Bogoliubov phonon, the tunneling properties of 
each component J„ and are more crucial than the quasiparticle tunneling current Jqp 
given by the sum of them. In Sees. Ill and IV, we show that each and has the 
same tunneling properties as those of supercurrent in the low energy region, leading to the 
anomalous tunneling of Bogoliubov phonon. 



III. SUPERCURRENT BEHAVIOR OF LOW-ENERGY BOGOLIUBOV MODE 

In this section, we first consider the tunneling of low-energy Bogoliubov phonon in the 
absence of supercurrent. In this case, the condensate wavefunction in Eq. (12.61) reduces to 



^'g=o(a^) = 7(3;) = tanh(|x| + Xq). 



(3.1) 



Since we are interested in the anomalous tunneling effect, we consider the low energy region 
{E <^1). In this regime, the Bogoliubov excitation spectrum has the linear dispersion E = p, 
where p (> 0) is the momentum of incident wave coming from x = —00. For the momenta 
of the four particular solutions in Eq. (12.171) . we may take (^1,^2,^3,^4) = {p, — 2i, —2i) 
within the accuracy of 0{p). Determining the coefficients {A, B,C, D) in Eqs. (I2.19P and 
(I2.20p . we obtain [2I, to the accuracy of 0{p), 



u<[x) 



Ur[X 



uAx) , 

vi{x) I \ Vr[X 




-ipx 



+ 



+ Pp 




(3.2) 





The coefficients (a, rj) are given by 



( 






a = 












< (3 = 


1 






il 


2 




V = 


a 





l-7o) 



270 



2 + 70 + 7o 
(l + 7o)7o ' 



l + 7o^ 



2 ' 



(3.3) 



(3.4) 
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where we have simply written 70 = 7(x = 0). In Eqs. (13.21) and (13. 3p . the incident wave 
{ui,vi) and transmitted wave {ut,vt) are obtained from the particular solution {ui,vi) in 
Eq. (I2.15p . while the reflected wave {ur,vr) is obtained from {u2,V2)- In the low energy 
limit, they are given by 




y 7(3;) - ^[1 - 7(3;)^] - - i] ^ 



(3.5) 





^^{x)+'-[l-^{xf]-^[^{x) + t] ^ 



(3.6) 



The localized component {ul,vl) in Eqs. (13. 2p and (13.30 is obtained from (u^jV^) and 
(M4,f4), which has the form 







["■■I 





-[l-7(x)']+p[l-7(x)] 
[l-l{xy]+p[l--f{x)] 



(3.7) 



We note that, in the last terms of Eqs. (13. 2p and (13.31) . the exponential factors e^^^ appearing 
in (^3,^3) and (M4, f4) have been absorbed into {ul,vl) by using the identity 



-2(|x[+xo)[i + ^(3.)] = l-7(x). 



(3.8) 



To show the supercurrent behavior of low-energy Bogoliubov phonon in the tunneling 
process, we expand Eqs. (13. 2p and (13.31) in terms of p to the first order. For example, M>(a:) 
then becomes 



f \ ( \ < ■ ^(^) - 70 , . . ^ 

Uy{x) = 'y[x) + ip \-ipx'y[x) 



+ ip 



l + 7o^ 



7o 

2li^) + 



70 . , P 70 ri .,/^\2i , P.J^^ , P. 



21+7o' 



Using the identity 



d_ 
dx 



7(x) = sgn(x) ^1 — 7(0;) 



we find that Eq. (13. 9p is equivalent to the expression. 



(x) = ^l+^e^^oe'P- l{^l+^x + T^TT^r^) + W 

L 2 i + 7n 



-70 
7o 



(3.9) 



(3.10) 



(3.11) 
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within the accuracy of 0{p). Expanding the condensate wave function in Eq. fl2.6p in the 
long wave length limit, we obtain (x > 0) 

^^(a;) = e'"" (j{x) + tq'^M^) . (3.12) 
^ 7o ^ 

Comparing Eq. (13. lip with (13.120 . we find that Uy{x) essentially has the same form as the 
condensate wavefunction \E'g=p(x) in the supercurrent state with momentum p. The same 



conclusion is also obtained for u<(x). Introducing = y/1 + px and 

1 -^^'T^^ 



u{y^) = —=e '+-^u{x = y+/^l+p), (3.13) 
VI +p 

we obtain 

%+) = §p(y+,y+). (3.14) 

Here, is given by Eq. (12. 6p . where xo in 7 is replaced by 

n = xo + ^^^p. (3.15) 

The same analysis is applicable to v{x). From Eqs. (I3.2p -( l3r7j) . Vp{x) is found to be 
related to Up{x) as Vp{x) = u_p{x)*. Using this, we obtain 

v{y^) = ^p{y^,Y^), (3.16) 

where = ^/T^^x, 



2(1 + li) 

and 



Y_ = xo- (3.17) 



70 



v{y.) = ^^=e '"^ovix = y.s/l^). (3.18) 

Equations (I3.14p and (I3.16P clearly show that, in the low energy region, u{x) and v{x) 
have the same properties as the condensate wavefunction in the supercurrent state. Namely, 
the currents J„ and J„ tunnel through the barrier without refiection in the low energy region. 
The tunneling of low energy Bogoliubov phonon described by the wavefunction (m, v) is also 
not accompanied by refiection, which 

The coefficient of the refiected wave component given by the second term in the right hand 
side of Eq. (13. 2p is proportional to so that the refiection probability is proportional to p^. 
Namely, deviation from the perfect transmission starts from O(p^). 
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Since u{x) and v{x) have the same form as the condensate wavefunction, the currents Ju 
and Jy are immediately obtained from Eq. (12.101) as 

Ju = -Jv = p. (3.19) 

Namely, and Jy are conserved everywhere in the low energy region. When we evaluate 
Jy and Jy of the incident wave, one finds Jy = — Jy = p + OijP') for a; — oo. This also 
confirms the perfect transmission of Jy and Jy, as well as Bogoliubov phonon. 

We note that, although Eq. (13.191) looks different from Eq. (I2.23P (where Jy and Jy are 
constant in the low energy limit), the reason is simply due to the assumed magnitude of the 
incident wave in Eq. (13. 2p . When we choose the magnitude of the incident wave (u/,f/) 
far away from the barrier so as to be equal to the plane wave solution (wp, Vp) given by 
Eq. (12.221) . Eq. (I3.19P is replaced by Jy = — Jy = 1/2. Under this normalization, a finite 
quasiparticle current Jqp obtained from higher order terms in J„ and Jy with respect t o p 



is proportional to p, as in the uniform system discussed in Sec. II. It has been shown |10| 
that this finite Jq is not conserved near the barrier due to the induction of supercurrent 
counterfiow. 

Figure [2] shows the incident, reflected, transmitted, and localized components of u{x) 
in Eqs. (13.21) and (13.31) . Although the perfect transmission looks as if the potential bar- 
rier is transparent for low-energy Bogoliubov phonon, the reflected component Ur{x) = 
iapuR{x)e~^'P^ , as well as the localized component Ul{x) = sgn(x)(2/2 + Pp)ul{x), actually 
contribute to the solution. Indeed, these components are necessary to satisfy the boundary 
conditions at a; = 0. However, since the coefficient of the reflected component Ur{x) is 
proportional to p, the plane wave factor e~*^^ in this component is actually irrelevant in 
the present treatment within 0{p). Thus, although Eq. (13.131) involves the reflected wave 
component, the reflected "plane wave" is not actually involved in it. 

We note that, while the localized state described by the last terms in Eqs. (13. 2p and 
(13. 3p does not carry current by itself, it still contributes to the tunneling current through 
the coupling with the propagating wave in the wavefunction. To see this, we divide Jy into 
contributions coming from each wave component and their couplings. Then we obtain 

Jy = Jp + JpL. (3.20) 

Here, Jp = lm[Up{x)*dxUp{x)] is the contribution from the right-going wave Up{x) = 
9{-x)ui{x)e'P'' + e{x){l - irip)uT{x)e'P'' , while Jpl = lm[Up{x)*d^ULix) + UL{xyd^Up{x)] 
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FIG. 2: Spatial variation of u(x) given by Eqs. (|3.2p and ()3.3p . We also show the incident 
wave component Ui{x) = ui{x)e^^^ , reflected wave Uji{x) = iapufi{x)e^^P''^ , transmitted wave 
Ut{x) = (1 — irip)uT{x)e'^P^ , and localized component Ul{x) = sgn(x)(i/2 + I3p)ul{x). We take 
p = 0.05 and xq = 0.5. 

describes coupling effects between the right-going wave and the localized state. The other 
current components involving the reflected wave Ur{x), such as Jr = lm[UR{x)*dxUji{x)], 
Jrp = lm[UR{x)*d^Up{x) + Up{x)*d^UR{x)], and Jrl = lm[t/jj(x)*(9^[/L(x) + 
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FIG. 3: Effects of localized state Ul{x) on J^. We take p = 0.05 and xq = 0.5. 
Ur^xYOxUl^x)]^ can be ignored in the low energy limit. Calculating Jp and Jpi, we obtain 



Figure [3] shows the contribution of Jp and Jpi to the current J„. When we ignore 
effects of localized state Ul{x), the current Ju (= Jp) is not conserved near the potential 
barrier, where the suppression of the condensate wavefunction \l/q=o(x) is remarkable. The 
enhancement of Jp is cancelled out by the counterfiow Jpi (< 0) originating from the 
coupling of the propagating wave Up{x) with the localized state Ul{x). 

Since u{x) and v{x) have the same form as the condensate wavefunction we can 

expect that the Bogoliubov equations (12.111) and (12.121) are also related to GP equation (12. Sp 
in the low energy region. To see this, it is convenient to note that the Bogoliubov equations 
(12. lip and (I2.12p can be formally written as 




(3.21) 




(3.22) 




(3.23) 



Here Wu and are given by 




(3.24) 
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In the absence of the supercurrent (g = 0), setting /i = 1, one can rewrite Eqs. (13.221) and 
(I3.23P in the form 

l^^Siy,) + '^)u(y,)=0, (3,25) 



2dyl VTTp l+p 

.i-i?l-l + -i=%_) + 5i^)*(,_) = 0. (3.26) 

2 dyl — p 1 — p / 

When Wu{x) / (1 +p) and Wy{x) / (1 —p) coincide with the nonhnear term in the GP equation, 
Eqs.f l3:25|l and fl3:26|l reproduce the GP equation (Q- Indeed, substituting Eqs. ([31])- ([31]) 
into Eqs. (13.241) . we obtain Wu{x)/{1 +p) = |£t(y+)p and W„/(l - p) = |-u(?/_)p, within the 
accuracy of 0{p). Thus, as expected, Eqs. (13.251) and fl3.26p have the same form as the GP 



17|. 



equation in the supercurrent state with momentum p 

We note that the potential barrier is modified as V/ y/1 ± p in Eqs. (I3.25P and (13.260 . 
which can explain the reason for the shift of xq given by Y± in Eqs. (I3.15P and (I3.17p . For 
example, when one solves the "GP" equation (I3.25p . the equation for is obtained as 

1 -7(0,F+)' = -^7(0, r+) ^ (l - |)V7(0,F+), (3.27) 

where we have ignored higher order terms than 0{p). In Eq. (13. 270 . we are using the 
notation 7(0, to emphasize that Xq is different from Y^. Noting that Xq is obtained from 
the equation 1 — 7(0, XqY = ^7(0, Xq), one can solve Eq. (13.271) within the accuracy of 0{p). 
The result is 

-^'"■-° + 2(i + °io.l)^) ^>- (^•^«) 

Here, we have used the identity in Eq. (I3.10p . Equation (13.28P reproduces F+ in Eq. (13.150 . 
In the same manner, one can confirm that the shift Yl of xq in v{y-) is due to the modified 
potential barrier V/ — p in the "GP" equation (13.260 . 



IV. QUASIPARTICLE TUNNELING IN THE SUPERCURRENT STATE 

In this section, we examine effects of supercurrent on the tunneling of Bogoliubov phonon. 
Assuming that the momentum q carried by condensate, as well as the incident momentum p 
of the Bogoliubov mode, are very small, we treat them within the first order. We also assume 
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that the supercurrent Jg is much smaller than the critical current. The critical supercurrent 
state is considered in Sec. V. 

Within the accuracy of 0(g), we may still ignore q in Eqs. f l2.7p and fl2.9p . Namely, 
7(x) and Xq are not modified by supercurrent. In addition, we can also set E = p and 
(Pi)P2,P3,P4) = (p, — p, 2i, — 2i) under the assumption of small p and q. The Bogoliubov 
wavefunction has the form 

0{-x)+ ^ 





[x)\ 






» ) 





where 





e*^^ + tap 




-vpx 



(4.1) 




(4.2) 





= (1 - irip) 1 


( ut{x) \ 


\ v>{x) ) 




[ vt{x) I 



^vpx 




(4.3) 



Here, coefficients {a,P,ri) and the localized component {ul,ul) are not affected by super- 
current, and they are given by Eqs. (13.41) and (13.71) . respectively. In contrast, the incident 
wave {ui,vi), reflected wave (m_r, fi?), and transmitted wave {ut,vt) are modifled to be 




^ li^) + ^[1 - 7(x)^] + |[7(x) + + lixf] ^ 
y 7(x) - ^[1 - lixf] - |[7(x) + lix)'] J 



(4.4) 




^ 7(^) - ^[1 - 7(^)'] + ?[7(a:) + + lixf] ^ 



2 

7(a;) + ^[l 



lyxj 



I [7(3;) 



.|[l + 7(x)^] 



(4.5) 




^ 7(x) - ^[1 - 7(x)^] + |[7(x) + lixf] ^ 

[ 7(^) + ^[1 - lixf] - ^bix) + + + lixf] J 



(4.6) 



To see the supercurrent behavior of low-energy Bogoliubov phonon in the supercurrent 
state, we expand Eq. (14. ip in terms of p and q to the flrst order. In this procedure, we also 
expand the phase factor e^*^' deflned in Eq. (12.80 as 



IT^- + 0{q'). 
7o 



(4.7) 
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For example, u{x > 0) becomes 

u{x > 0) = 7(0;) + i{p + q) h i{p + q)x'^{x) 

7o 



Comparing Eq. (14. 8p with Eq. (13.91) . we find that Eq. (14.81) can be written as, within the 
accuracy of 0(p) and 0(g), 



u{x > 0) = yrT^e'^^e*(P+^)- 7(^1+^0; + ,J + i{p + g) ^^"^ . (4.9) 

L 2(1 + 7o) 7o J 

This is essentially the same form as the condensate wavefunction \l/p+g(x) in (I3.12p . We also 
reach the same conclusion for m(< 0). As a result, we obtain 

u{y+) = ^,+,{y+,Y+). (4.10) 

Using the relation Vp{x) = u*_p{x), we also find 

^(y_) = ^p_,(i/_,r_). (4.11) 

Equations (I4.10p and (14. lip indicate that tunneling properties of the current J„ and 
are the same as those of supercurrent in the low energy region. Namely, the Bogoliubov 
phonon still shows the perfect transmission in the supercurrent state [4]. Equations (I4.10p 
and (14. lip also show that the supercurrent Jg = q affects J„ and J^, as J„ = p + g and 
Jy = —p + g. This implies that, in the critical supercurrent state, the "supercurrent" J„ or 

exceeds the critical current, leading to the breakdown of the perfect transmission. We 
will confirm this in Sec. V. 

We note that Eqs. (I3.25P and (I3.26P are also valid for the supercurrent state. Substituting 
Eqs. fICTD and KT\\ into Eq. ([321, we obtain PV„/(1 + p) = \u{y+)\^ and Wj{l~p) = 
|{;(?/_)p. Thus, Eqs. (I3.25P and (13.260 reduce to the GP equation, as expected. 

V. BREAKDOWN OF ANOMALOUS TUNNELING EFFECT IN THE CRITICAL 
SUPERCURRENT STATE 

In Ref. 4], the breakdown of the anomalous tunneling is predicted in the critical super- 
current state, when the barrier V is high enough so that the tunneling of supercurrent can 
be regarded as the Josephson effect. In this section, we also consider the same situation 

16 



to see if the absence of the anomalous tunnehng in the critical supercurrent state can be 
understood as that the currents Ju or exceeds the critical value. 

We briefly summarize the critical supercurrent Jc (= (jc) in the Josephson coupling regime. 
To derive the current-phase relation in this regime, it is convenient to write the condensate 
wavefunction '^q{x) in the form 

^,{x) = v/g2 + 7(x)V'?"e'^s°(")^. (5.1) 

Here, the phase (f){x) is defined by 

= 2ftan-i^ -tan-i^y (5.2) 
V q q J 

Equation (15.21) can be rewritten as 



$(x) 

g = 7otan^— , (5.3) 

where $(x) = (pix) — 2[tan~-'^(7(a;)/g) — 7r/2]. When the barrier potential V is very large, q 
and 7o are very small and one can expand them in terms of J]. Within the accuracy of 
0{V-^), Eq. ([23]) gives 

70 = ^70+?')- (5.4) 

From Eqs. (15. 3p and (15.41) . one finds 

o = sin <1>, 

^ 2^ (5.5) 
70 = 2^7(1 + cos $). 

Since the supercurrent is uniform, the x dependence of $(a;) in Eq. (15.51) is actually irrele- 
vant. When we take the phase 0(x) at a; = 00 (= (po), we obtain the well-known Josephson's 
current-phase relation, Jg = {1/2V) sin^o within the accuracy of 0(y^^)[^]. Equation (15. 5p 
gives the critical current Jc in the Josephson coupling regime as 

Jc = qc = (5.6) 

Now we construct the Bogoliubov wavefunction {u,v). In the following discussion, we 
always assume p ^ g ~ 70 ~ 1/2V^ <^ 1. The Bogoliubov wavefunction, satisfying the 
boundary conditions, are given by Eqs. (I4.ip - (I4.6I) . where the coefficient (3 in Eq. (14. 3 p is 
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replaced by [3'. The coefficients {a,P,P',ri) are given by 

1 .70 



"=2 



V 



1 1 - 2g 70 



4 



, , 7o-9'' 
1 _ 1 - g 7o 

"4 2 jl-q^' 
--(1-4.)-^ 



(5.7) 



2 • 

We note that Eq. fl5.5p gives 7o = Q'c = l/2\^ in the critical supercurrent state. In Sec. IV, 
we have expanded the Bogoliubov wavefunction {u, v) in terms of p and q, where we have 
imphcitly assumed g/70 <^ 1. (See, for example, Eq. 04.71) .) In the present case (J^ < Jc), 
although we can still assume g <^ 1, we have to treat 70 as the same order as q. Namely, 
the expansion in terms of q/'~fo, such as Eq. (14.71) . is not allowed in the present case. In 
addition, because of g ~ 70, we have to carefully treat the factor 1/(79 ~ 0'^) ^'i- (15-71) . 
However, even in this case, we can show that u and v still have the supercurrent properties 
near the barrier at x = 0. 

To show this, we ffist consider u{x). Within the accuracy of 0{x), u{x > 0) reduces to 

7o 



u{x > 0) = ( 1 + ip „ ) ( 7(x) - iQ+ - 2pq- 



,2' 



7o - 

where = p + q. We also expand the factor e"*^'' in Eq. (14.11) to 0{p), which gives 

lo + iq 7o + zQ+ . 7o 



-tea 



Using Eqs. ([EHD and ([ESD, we obtainQ 



1 — ip 



ll + q' 



(5.8) 



(5.9) 



u{x > 0) 



7(o,^+) + 2g+ 



V7(0, Z+)2 + Ql 
where 7(x, Z+) equals 7(x) where Xq is replaced by 

7(0, Z. 



e^«+^[7(x,Z+)-zgH 



(5.10) 



Z+ = Xo- 2pq- 

7(0, Z 

In the same way, we can rewrite ^(x < 0) as 

7(o,^-)-^g+ 



u{x < 0) 



e*«+^[7(x,Z+)+zQ. 



v/7(o, z+y + Q\ 

Comparing Eqs. (15.101) and (15.121) with Eq. (12. 6p we find near the barrier (x <^ 1) 

U{X) = ^p + g(X,Z+). 



(5.11) 



(5.12) 



(5.13) 
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Applying the same analysis to v{x), we obtain 

vix) =^p.,ix,Z_), (5.14) 

where 

Z_ = Xo + 2pq (5.15) 
7(0, Z_) - 

Thus, we conclude that u and v still have the same properties as those of the condensate 
wavefunction near the barrier at a; = 0. Their currents are given by = p + q and 
= —p + g at X <^ 1. 

Using the fact that 70 satisfies Eq. (15. 4p . we obtains, within the accuracy of 0{p) 

7(0,Z±) = l^(7(0,Z±)2 + Q|). (5.16) 
We define the phase (p±{x) (a; -C 1) by 

^.W = 2(ta„-I<|l^-tan-'2(MJ). (5.17) 

This is analogous to the phase 0(x) introduced in Eq. (15.21) . From Eqs. (I5.16P and (15.171) . 
we obtain 

Q± = ^ sin$±, 

(5.18) 

7(0, Z±) = ^(l + cos$±), 

where $±(x) = 4>±{x) — 2[tan^"'^(7(x, Z±)/Q±) — 7r/2]. Since the currents J„ = = p + q 
and Ju = —Q- = —p + q are uniform near the barrier, the x dependence of $± in Eq. (I5.18P 
can be actually ignored as far as x ^ 1. Equation (15.180 shows that both and Jy must 
be smaller than the upper limit of the supercurrent Jc = 1/2V. However, this condition 
cannot be satisfied in the critical supercurrent state {q = Jc = ^/2V). Namely, J„ or J„ 
always exceeds Jc, depending on the direction of the momentum p of Bogoliubov phonon. 
As a result, the perfect transmission is not obtained in the critical supercurrent state. We 



breakdown of the anomalous tunneling 



emphasize that this gives simple explanation for the 
in the critical supercurrent state predicted in Ref. {4]. 

We note that the coefficients {a, /?, 7) in Eqs. (15. 7p diverge when g = 70 = 1/2V^. This 
means that the components u and v in the Bogoliubov wavefunction no longer have the 
same form as the condensate wavefunction in the critical supercurrent state. In this case, 
u and V would have different forms from the condensate wavefunction, leading to a finite 
refiection probability at the barrier. 
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VI. SUMMARY 



In this paper, we have investigated tunnehng properties of low-energy Bogohubov ex- 
citations in a Bose superfluid. Using the exactly solvable tunneling problem through a 
5-functional potential barrier, we have clarified that, in the low energy region, the com- 
ponents u{x) and v{x) in the Bogoliubov wavefunction {u{x),v{x)) have the same form as 
the condensate wavefunction in the supercurrent state. Based on this result, we have given 
physical explanation for the anomalous tunneling effect of Bogoliubov phonon, as well as 
the absence of this phenomenon in the critical supercurrent state, in a consistent manner. 

The current and Jy associated with u{x) and v{x), respectively, tunnel through the 
barrier without reflection, as in the case of supercurrent tunneling. This gives physical 
explanation for the perfect transmission of the Bogoliubov phonon predicted recently. We 
also showed that the supercurrent behaviors of J„ and J^, as well as the resulting perfect 
transmission of Bogoliubov phonon, are also obtained in the supercurrent state unless the 
magnitude of the supercurrent reaches the critical value. 

In the Josephson coupling regime (which is realized when the tunneling barrier is very 
high) , the upper limit of and Jy is given by the critical current Jc of the ordinary super- 
current. In the critical supercurrent state, Ju or always exceeds Jc, irrespective of the 
direction of the momentum of Bogoliubov phonon. Because of this, the perfect transmission 
of Bogoliubov phonon is no longer obtained in the critical supercurrent state. This can ex- 
plain the reason why the breakdown of the anomalous tunneling effect occurs in the critical 
supercurrent state. 

In this paper, we have used a simple 5-functional potential barrier. While this model 
enables us to treat the GP and Bogoliubov equations analytically, any realistic potential 
barrier should have a finite potential width. However, although our model is simple, we still 
expect that the essence of our results in this paper would be valid for the case with a more 
realistic barrier potential. We will extend our study to more general and realistic cases in 
our future paper. Since the Bogoliubov phonon is one of the most fundamental phenomena 
in the superfiuid phase, the observation of the anomalous tunneling would be important in 
understanding basic properties of this collective mode. 
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